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We examine statistical isotropy of large scale anisotropies of the Internal Linear Combination 
(ILC) map, based on three year WMAP data. Our study reveals no significant deviation from 
statistical isotropy on large angular scales of 3-year ILC map. Comparing statistical isotropy of 
3-year ILC map and 1-year ILC map, we find a significant improvement in 3-year ILC map which 
can be due to the gain model, improved ILC map processing and foreground minimization. 



I. INTRODUCTION 

The Cosmic Microwave Background (CMB) anisotropy 
has been shown to be a very powerful observational probe 
of cosmology. Detailed measurements of the anisotropies 
in the CMB can provide a wealth of information about 
the global properties, constituents and history of the 
Universe. In standard cosmology, the CMB anisotropy 
is expected to be statistically isotropic, i.e., statistical 
expectation values of the temperature fluctuations (and 
in particular the angular correlation function) are pre- 
served under rotations of the sky. This property of CMB 
anisotropy has been under scrutinyafter the release of 
the firstyear of WMAP data Q IlJiJ3, IES EE 

IH El El III El 13 III 13 IlllflMlllMllir We 
use a method based on bipolar expansion of the two point 
correlation function which is an improved and enhanced 
follow up on our previous work on first year WMAP data 
pil |25| . This method is shown to be sensitive to struc- 
tures and patterns in the underlying two-point correla- 
tion function. 

We apply our method to the improved Internal Linear 
Combination (ILC) map 26] , based on three year WMAP 
data[5Cj. We choose the ILC map for testing statistical 
isotropy (SI) of the CMB anisotropy for the following 
reasons 

1. The ILC is a full-sky map and hence is easier 
to work with. Masking the sky results in viola- 
tion of statistical isotropy. An originally SI CMB 
anisotropy map deviates from SI after masking |25j| . 

2. Residuals from Galactic removal errors in the three- 
year ILC map are estimated to be less than 5 [iK 
on angular scales greater than ~ 10 deg |26|. Hence 
at low-Z, multipoles are not significantly affected by 
foregrounds. In addition, it is interesting to exam- 
ine the above statement by testing the statistical 
isotropy of the ILC map. 

3. On large scales, the three- year ILC map is believed 
to provide a reliable estimate of the CMB signal, 



with negligible instrument noise, over the full sky 

m 

These properties of the ILC map allow us to study the 
cosmological signal on large scales. In addition, there 
are theoretical motivations for hunting for SI violation 
on large scales of CMB anisotropy. Topologically com- 
pact spaces lia. 1471 liM |49l and anisotropic cosmological 
models [30t l3ll 13211331 Im , I3M l3r| are examples of this. 
Both observational artifacts and the above theoretical 
models cause a departure from statistical isotropy and it 
has been shown that our method is a useful tool to find 
out these deviations (see e.g. [3E1I13)- The rest of 
this paper is organized as follows: Section [D] is a brief 
introduction to temperature anisotropy of CMB. Section 
^Jdescribes the formulation of statistical isotropy in gen- 
eral. Section HVI is a description of estimators we use to 
look for deviations from statistical isotropy. We present 
the application of our method on the WMAP data in 
Section [V] Section IVII contains discussion on the cos- 
mological implications of our null detection of deviations 
from statistical isotropy on large angular scales in 3-year 
ILC map of WMAP data, and in Section lYlIl wc summa- 
rize our results. 



II. CHARACTERIZATION OF CMB 
TEMPERATURE ANISOTROPY 

The CMB anisotropy is fully described by its tem- 
perature anisotropy and polarization. The temperature 
anisotropy is a scalar random field, AT(n) = T(h) — T , 
on a 2-dimensional surface of a sphere (the sky), where 
n = (9, <j)) is a unit vector on the sphere and To — 
J ^ L T(n) represents the mean temperature of the CMB. 
It is convenient to expand the temperature anisotropy 
field into spherical harmonics, the orthonormal basis on 
the sphere, as 



AT(n) 



E 

l,m 



aimYlrnifl) 



(1) 
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where the complex quantities, ai m are given by 



ai m = / da h Y* m (n)AT(n). 



(2) 



Statistical properties of this field can be characterized by 
n-point correlation functions 



(AT(ni)AT(n 2 ) • • • AT(h n )). 



(3) 



Here the bracket denotes the ensemble average, i.e. an 
average over all possible configurations of the field. CMB 
anisotropy is believed to be Gaussian [ifl Elj . Hence the 
connected part of n-point functions disappears for n > 
2. Non-zero (even-n)-point correlation functions can be 
expressed in terms of the 2-point correlation function. As 
a result, a Gaussian distribution is completely described 
by the two-point correlation function 



C(n,ri) = (AT(n)AT(n')). 



(4) 



Equivalently; as it is seen from eqn. for a Gaussian 
CMB anisotropy, ai m are complex Gaussian random vari- 
ables too. Therefore, the covariance matrix, (a/ m a*/ m /), 
fully describes the whole field. Throughout this paper we 
assume Gaussianity to be valid. 



III. STATISTICAL ISOTROPY 

Two point correlations of CMB anisotropy, C(hi, 712), 
are two point functions on S 2 x S 2 , and hence can be 
expanded as 

C(ni, n 2 ) = A iM\hi 2 Y iM 2 (ni, ri 2 ). (5) 



Here A 



iMlhk 



are coefficients of the expansion (hereafter 

BipoSH coefficients) and Y^ 2 (fii, n<z) are bipolar spher- 
ical harmonics defined by eqn. i|AI|l . Bipolar spherical 
harmonics form an orthonormal basis on S 2 x S 2 and 
transform in the same manner as the spherical harmonic 
function with I, M with respect to rotations |42j. One 
can inverse-transform C(ni, n?) in eqn. © to get the 
coefficients of expansion, A eM n l2 , ky multiplying both 

sides of eqn.® by Y£, M ? (ni , ^2) and integrating over all 
angles. Then the orthonormality of bipolar harmonics, 
eqn. 1A2II . implies that 



A £M\hi 2 = J dSlnx J dttf l2 C(h 1 ,h 2 )Y*l} l ' 2 (h 1 ,h 2 ). 

(6) 

The above expression and the fact that C(h\, 712) is sym- 
metric under the exchange of n\ and h 2 leads to the fol- 
lowing symmetries of A lM ^^ 



A 



eM\hh 



(_i\(h+l2-L) a 



(J) 



A IM\U — A IM\U ^,2fe, 



k = 0, 1, 2, 3, 



It has been shown |25( that Bipolar Spherical Harmonic 
(BipoSH) coefficients, ^gM\hi 2 ^ are * n ^ ac * l mear combi- 
nations of off-diagonal elements of the covariance matrix, 
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(8) 



where C, „ , „ are Clebsch-Gordan coefficients (see 
the Appendix). This clearly shows that A £M ^ il2 com- 
pletely represent the information of the covariance ma- 
trix. When statistical isotropy holds, it is guaranteed 
that the covariance matrix is diagonal, 



{ai m a v 



Cl Su'S„ 



(9) 



and hence the angular power spectra carry all informa- 
tion of the field. Substituting this into eqn. JSJ) gives 



A 



CM\IV 
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MO- 



(10) 



The above expression tells us that when statistical 
isotropy holds, all BipoSH coefficients, A 



£M\W 



are zero 



except those with I = 0, M = which are equal to the 
angular power spectra up to a (— l) 1 (21 + 1) 1 / 2 factor. Bi- 
poSH expansion is the most general way of studying two 
point correlation functions of CMB anisotropy. The well 
known angular power spectrum, Ci is in fact a subset of 
the corresponding BipoSH coefficients, 



a 



(-1)' 



--A, 
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(11) 



Therefore to test a CMB map for statistical isotropy, one 
should compute the BipoSH coefficients for the maps and 
look for nonzero BipoSH coefficients. Statistically signif- 
icant deviations from zero would mean violation of sta- 
tistical isotropy. 



IV. ESTIMATORS 

Given a CMB anisotropy map, one can measure Bi- 
poSH coefficients bye the following estimator[5lj 



.4 



IM\IV 



IM 

Irnl' m' 



(12) 



where W\ is the Legendre transform of the window func- 
tion. The above estimator is a combination of C; and 
hence is un-biased|52|. However it is impossible to mea- 
sure all ^4 M /|;;/ individually because of cosmic variance. 
Combining BipoSH coefficients helps to reduce the cos- 
mic variance |53J. There are several ways of combining 
BipoSH coefficients. Here we choose two methods. 
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which statistical isotropy is valid, BipoSH coefficients are 
given by eqn and therefore SI condition implies a 

null BiPS, i.e. Kg = for every t > 0, 



K 5. 



to ■ 



(14) 



Non-zero components of BiPS imply break down of sta- 
tistical isotropy, and this introduces BiPS as a measure 
of statistical isotropy. It is worth noting that although 
BiPS is quartic in ai m , it is designed to detect SI vio- 
lation and not non-Gaussianity 0, |2^| Q . An 
un-biased estimator of BiPS is given by 



E 

IV M 



A 



.'A/I//' 



58, 



(15) 



where 25^ is the bias that arises from the SI part of the 
map and is given by the angular power spectrum, Ci, 



23/ 



(hf) 
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(16) 



t+h 
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The above expression for 58^ is obtained by assuming 
Gaussian statistics of the temperature fluctuations |25j 

0. 



B. New Method: Reduced Bipolar Coefficients 

The BipoSH coefficients of eqn. I|12|) can be summed 
over I and I' to reduce the cosmic variance, 



FIG. 1: Top: A bipolar map generated from bipolar coeffi- 
cients, Aim, of 3-year ILC map. Middle: bipolar map based 
on 1-year ILC map and Bottom: differences between the 
two maps (note the scales). The top map (ILC-3) has smaller 
fluctuations comparing to the middle one (ILC-1) except for 
the hot spot near the equator. Differences between these two 
maps mostly arise from a band around the equator in bipolar 
space. Both ILC maps are smoothed by a band pass filter, 
W s (k =2,l s = 10). 



A. First Method: BiPS 

Among the several possible combinations of BipoSH 
coefficients, the Bipolar Power Spectrum (BiPS) was 
proved to be a useful tool with interesting features |43| . 
BiPS of CMB anisotropy is defined as a convenient con- 
traction of the BipoSH coefficients 



K( 



E 

1,1', M 



A 



IM\W 



> 0. 



(13) 



BiPS is interesting because it is orientation independent, 
i.e. invariant under rotations of the sky. For models in 



.4 



IM 
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(17) 



1=0 l' 



These reduced bipolar coefficients, Aim, by definition re- 
spect the following symmetry: 



A 



IM 



{-1) M A} 



-Mi 



(18) 



which indicates Aio are always real. When SI condition 
is valid, the ensemble average of Aim vanishes for all I 
and M 



{A 



(Ml 



0. 



(19) 



In any given CMB anisotropy map, Aim would fluctuate 
about zero. A severe breakdown of statistical isotropy 
will result in huge deviations from zero. Reduced bipo- 
lar coefficients are not rotationally invariant, hence they 
assign direction to the correlation patterns of a map. We 
can combine Aim further to define a power spectrum 
similar to how ai rn are combined to construct the angu- 
lar power spectrum, C;. We define 



1 



2^+1 



E A tMA iM . 



(20) 
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FIG. 2: Real part of Aim's of ILC-3 (red square points) and 
ILC-1 (blue stars) for a W S (10,2) filter that roughly keeps 
multipoles between 2 and 15. I and M indices are combined 
to a single index n = £(£+l)+ M + 1 and the blue dotted lines 
define 1-a error bars derived from 1000 simulations of SI CMB 
anisotropy maps. Almost all of Aim's of ILC-3 are smaller 
than those of ILC-1 which means ILC-3 is more consistent 
with statistical isotropy. 



The above estimator is rotationally invariant. It has a 
positive bias and hence it has similar issues that have 
been addressed for the BiPS studies earlier. This means 
although the ensemble average of Aim for a statistically 
isotropic case is zero, ensemble average of Di is always 
greater than zero. However a major deviation from sta- 
tistical isotropy will result in a big Di (compared to that 
of a SI case). In section [V] we compare Aim of ILC 
map against an average of 1000 simulations of statisti- 
cally isotropic maps. We defer detailed studies of D t to 
the future publication. 



V. APPLICATION TO THE WMAP DATA 

We carry out our analysis on 3-year ILC map and com- 
pare it to 1-year ILC map. In order to attribute a statis- 
tical significance to our results, we compare our results 
to 1000 simulations of SI CMB maps. a; m 's of these 
maps are generated up to an l max of 1024 (correspond- 
ing to HEALPixj^l resolution N sil i e = 512). Since we 
are only interested in large angular scales we smooth all 
maps with appropriate filters to cut the power on small 
angular scales. These filters are low pass Gaussian filters 



= N G exp 



21 + 1 
217+1 . 



(21) 



that cut power on scales (I > l s ) and band pass filters of 
the form 



W? = 2N S 



1- Jo 



2Z + 1 
2/t + l 
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FIG. 3: Fraction of Aim's in 1000 simulations that are smaller 
than Aim's of ILC. Several deviations in ILC-1 (red square 
points) have been corrected in ILC-3 (green impulses). Only 
deviations above 68% are shown. Almost all of Aim's of ILC- 
3 are smaller than those of ILC-1 which means ILC-3 is more 
consistent with statistical isotropy. Only real parts are shown. 
Points are always in pairs because of the symmetry in re- 



duced bipolar coefficients, Aim = (- 
smoothed with a W s (10,2) filter. 



l A*, 



All maps are 



that keep the power on scales corresponding to l t < I < 
l s -Jo is the bessel function and N G and N s are normal- 
ization constants chosen such that, 



y, (21 + l)Wi = ! 



I 



21(1 + 1) 



(23) 



i.e., unit rms for unit flat band angular power spectrum 

°' - 1(1+1) ■ 

for 3- 



We compute the BipoSH coefficients, A lM > u ,, 
year ILC map (ILC-3) for several window functions using 
eqn. I|12|). We combine these coefficients using eqn. <|17[) 
to obtain Aim- An interesting way of visualizing these 
coefficients is to make a map from them. Making a map 
from Aim is simply done similar to making a temperaure 
anisotropy map from a given set of spherical harmonic 
coefficients, aim) 



e (")=H A lM Y lM (h). 



(24) 



C=0 M=- 



(22) 



The symmetry of reduced bipolar coefficients, eqn. (T5J, 
guarantees reality of Q(n). The "bipolar" map based on 
bipolar coefficients of ILC-3 is shown on the top panel of 
Fig. ^ The map has small fluctuations except for a pair 
of hot and cold spots near the equator. To compare, we 
have also made a bipolar map of 1-year ILC map (ILC-1) 
from bipolar coefficients of ILC-1 (middle panel of Fig. 
[TJ. The difference map (Fig. ^ (bottom)) shows that 
differences between these two maps mostly arise from a 
band around the equator in bipolar space. As it is seen 
in Fig. 2] the bipolar map of ILC-3 has less fluctuations 
comparing to that of ILC-1. This is because almost all 
of A^m's of ILC-3 are smaller than those of ILC-1 (i.e. 
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are closer to zero). Reduced bipolar coefficients of the 
above maps are in Figure [21 in which £ and M indices 
are combined to a single index n — £(£ + 1) + M + 1 (only 
real part of Aim is plotted.). And the blue dotted lines 
define l-a error bars derived from f 000 simulations of SI 
CMB anisotropy maps. As it can be seen many spikes 
presented in Aim's of ILC-1 have either disappeared or 
reduced in ILC-3 (e.g. those around n = 20, 40 and a big 
spike at n = 111). To get a quantitative description of 
differences between ILC-3 and ILC-1 we compare them 
against 1000 simulations of SI CMB anisotropy maps. A 
simple x 2 comparison of Aim with simulations gives us 
a rough estimate of overall differences between the two 
ILC maps: ILC-3 has a smaller \ 2 than ILC-1. For a 
W S (10, 2) filter, the reduced \ 2 falls from 1.089 for ILC- 
1 to 0.9619 for ILC-3. Although \ 2 statistics is simple, 
it should be used with caution because it is only valid 
if every Aim is independent has a Gaussian distribution 
function. In order to study deviations of Aim from zero 
without worrying about the Gaussianity of the Aim, we 
look at the most deviant (biggest) A eM - We compare 
the biggest Aim's of ILC to Aim's of 1000 simulations to 
find out what fraction of simulations have Aim 's smaller 
than those of ILC maps. Figure shows the results. The 
horizontal axis is n = £(£ + 1) + M + 1 and the vertical 
axis is the fraction of Aim's in 1000 simulations that 
are smaller than Aim of ILC. In this figure red squares 
represent the ILC-1 while ILC-3 is represented by green 
lines. When a green line crosses a red point, Aim's of 
ILC-3 are greater than ILC-1, otherwise red points above 
green spikes show smaller Aim's for ILC-3. The results 
are interesting: several deviations in ILC-1 have been 
corrected in ILC-3. Specially on the largest scales, several 
deviations beyound 95% in ILC-1 have gone away in ILC- 
3 (red points above the blue dotted line in Figure [21 have 
been replaced by significantly smaller values). However 
there are a couple of exceptions that could be responsible 
for the hotter spot in bipolar map of ILC-3. 



Combining the BipoSH coefficients to construct bipo- 
lar power spectrum allows further examinations of ILC 
maps for departures from SI. We compute the BiPS us- 
ing eqn. I|15|) . It is worth mentioning that BiPS in this 
paper has been computed in a slightly different way than 
in our previous paper p4| . Here we compute the BiPS 
using eqn. I|15f) and we use the derived Ci from each map 
to estimate the bias, 25^, using eqn. Ijlb|l [55l|. The bias 
corrected BiPS is then averaged over 1000 simulations 
and is compared to bias corrected BiPS of ILC maps. 
BiPS results shown in Figure0]agree with our results on 
Aim- It can be seen that ILC-3 has a smaller bipolar 
power spectrum than ILC-1 and is more consistent with 
statistical isotropy. The same is true for Di estimator 
defined by eqn. (|20|) which we defer to the future publi- 
cations. We should emphasize that these results are only 
for large angular scales, Z < 25, and not beyond that. 
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FIG. 4: Bipolar power spectrum (BiPS) of the two ILC maps 
compared to average BiPS of 1000 simulations of statistically 
isotropic CMB maps. Both ILC maps are smoothed with 
a VF S (10, 2) window function which roughly retains multi- 
poles in the range of 2 < Z < 15. ILC-3 shows smaller BiPS 
than ILC-1, which means it is more consistent with statistical 
isotropy. Filtering the data with other functions show almost 
the same results. 



VI. DISCUSSION AND CONCLUSIONS 

The null results of search for departure from sta- 
tistical isotropy has implications for the observation 
and data analysis techniques used to create the CMB 
anisotropy maps. Observational artifacts such as non- 
circular beam, inhomogeneous noise correlation, residual 
striping patterns, and residuals from foregrounds are po- 
tential sources of SI breakdown. Our null results con- 
firm that these artifacts do not significantly contribute 
to large scale anisotropies of 3-year ILC map. We have 
also quantified the differences between 1-year and 3- 
year ILC maps. It is shown that 3-year ILC map is 
"cleaner" than 1-year ILC map at I < 25. This can be 
due to the gain model and improved ILC map process- 
ing and foreground minimization. It has also been ob- 
served that at large Z deviations from statistical isotropy 
occur which we think is because of residuals from fore- 
grounds. However we limit ourselves to the low-Z limit 
because in addition to observational artifacts, there are 
theoretical motivations for hunting for SI violation on 
large scales of CMB anisotropy. Topologically compact 
spaces Efl Wft EiM anc ^ anisotropic cosmological mod- 
els [30l I311 13a. I33L I34L l35l l36j are examples of this. Each 
of these models will cause departures from statistical 
isotropy in CMB anisotropy maps. And a null detec- 
tion of departure from statistical isotropy at low Z in the 
WMAP data can be used to put constraints on these 
models. Our measure is sensitive to axial asymmetries in 
the two point correlation of the temperature anisotropy 
[33. And this is even more significant now because the 
new measure of reduced bipolar coefficients does retain 
directional information. Our analysis doesn't show a sig- 
nificant detection of an "axis of evil" in the WMAP data. 
We have redone our analysis on ILC map filtered with a 
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low-pass filter that only keeps I = 2, 3,4 to search for a 
preferred direction at low multipoles. We have not been 
able to detect any significant deviation from statistical 
isotropy using various filters. We could not test the ef- 
fect of alignment of low multipoles on statistical isotropy 
because we had no theory or model to explain them. Va- 
lidity of statistical isotropy at large angular scales can put 
tight constraints on anisotropic mechanisms that are can- 
didates of explaining the low quadrupole of the WMAP 
and COBE data. It is worth noticing that our method 
can be extended to polarization maps of CMB anisotropy. 
Analysis of statistical isotropy of full-sky polarization 
maps of WMAP are currently under progress and will 
be reported in a separate publication. 

VII. SUMMARY 

We examine statistical isotropy of large scale 
anisotropics of the improved Internal Linear Combina- 
tion (ILC) map, based on three year WMAP data. In 
order to attribute a statistical significance to our results, 
we use 1000 simulations of statistically isotropic CMB 
maps. We have done our analysis using a series of fil- 
ters that span the low-^ multipoles. We only explicitly 
present the results for one of them that roughly retains 
power in the multipoles between 2 and 15. This reveals 
no significant deviation from statistical isotropy on large 
angular scales of 3-year ILC map. Comparing statistical 
isotropy of 3-year ILC map and 1-year ILC map, we find 
a significant improvement in 3-year ILC map which can 
be due to the gain model and improved ILC map pro- 
cessing and foreground minimization. We get consistent 
and similar results from other filters. 
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APPENDIX A: USEFUL MATHEMATICAL 
RELATIONS 



Bipolar spherical harmonics form an orthonormal basis 
of S 2 x S 2 and are defined as 



Ym(ni,n 2 ) = C hniihni2 Y h mi {ni)Y hm2 (h 2 ), 

(Al) 

in which Cf^ x TO2 are Clebsch-Gordan coefficients. 
Clebsch-Gordan coefficients are non-zero only if triangu- 
larity relation holds, {hfoi}, and M = mi + m 2 . Where 
the 3,7-symbol {abc} is defined by 



Jl if a + b + c is integer and \a — b\ < c < (a + b), 
\ otherwise, 



Orthonormality of bipolar spherical harmonics 
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